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Abstract 

We introduce a family of Jastrow pair product states for quasi one- 
dimensional spin systems. Depending on a parameter they interpolate be- 
tween the resonating valence bond ground state of the Haldane-Shastry model 
describing a spin liquid and the (dimerized) valence bond solid ground states 
of the Majumdar-Ghosh spin chain. These states are found to form an ex- 
cellent basis for variational studies of Heisenberg chains with next nearest 
neighbour interaction and bond alternation as realized in the spin-Peierls sys- 
tem CuGeOa- 
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I. INTRODUCTION 



Following the recent discovery of a spin-Peierls transition in the inorganic compound 
CuGeOj there has been growing theoretical interest in this instability of one-dimensional 
spin chains. It has been proposed in Refs. that frustrating next nearest neighbour 
(nnn) interactions in addition to an explicitely broken translational invariance due to lattice 
dimerization is necessary to obtain a consistent description of the experimental data. This 
leads to the following spin-| Hamiltonian 

N-l 

^ = E ((1 + (-l)'^)S, ■ S,+i + aS, ■ S,+2) . (1.1) 

j=0 

For 6 = 0, the model is invariant under translations by one lattice site. This case has been 
investigated in detail and includes the nearest neighbour Heisenberg chain, a = 0, where 
the complete spectrum can be obtained by means of the Bethe AnsatzS: the ground state 
is that of a spin liquid, has a vanishing spin gap and algebraically decaying correlations at 
T = 0. For 6 = 0, a = ^ the Hamiltonian (|1 . 1|) becomes that of the Majumdar-Ghosh (MG) 
modeli: here the system has a gapi and the exact ground state is known to be a product of 
nearest neighbour singlet pairs showing a twofold degeneracy 

« [0 ; 1] ■ ■ ■ [{N- 2) ; (N- 1)] , l^P^^l) (X [1 ; 2] ■ ■ . [{N- 1) ; 0] (1.2) 

{[a;b] = -^{\ t)a| l)b — I l)a\ T)f)) denotes the singlet state formed by the spins on sites a 
and b). For intermediate values of the nnn interaction the model (|1.1| ) remains gapless for 
a < Oc ~ 0.241lii'i with its low energy sector described by an effective level k = 1 SU(2) 
Wess-Zumino-Witten (WZW) conformal field theory. Increasing a beyond the 'conformal 
point' ac the nnn coupling becomes marginally relevant, producing an exponentially small 
gap A oc exp(const./(a; — ac))- In this phase the system is spontaneously dimerized. The 
properties of the system for a > 1 have recently been discussed in Ref. Extending the 
discussion to general couplings in the a-6 plane the system has gap above a dimerized 
ground state for any nonzero 5,0 with the valence bond ground states (|1.2| ) on the line 
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2a + 5 = lE. Going to larger S the Hamiltonian ( p..!]) corresponds to a ladder system of two 
coupled Heisenberg chain^. In addition to numerical studies, several mean-field theories 
have been proposed for the Hamiltonian ( |l.l| )Ei. In this framework physical properties 
have been calculated giving reasonable agreement with the experimental data available for 
energy gaps, Raman spectra and the susceptibility of CuGeOa and hence further support 
the Hamiltonian (|1 . 1|) as a model for this substance. 

In this paper we propose a family of variational states for the ground state and low lying 
triplet of the Hamiltonian (|1 . 1| ) . While it will not be possible within this variational approach 
to compute the thermodynamical quantities mentioned above, it turns out that the states 
proposed are excellent approximations to the true ground state of the system throughout 
the parameter region of interest here, namely 2a + 6 < 1. Hence they can be used to obtain 
very good variational bounds on energies and may provide a better understanding of the 
role of quantum fluctuations in these systems. 

In fact, variational states have already been applied successfully to some of the systems 
mentioned above: the Jastrow pair product wave function 

M M 

1^0^'^'^)= E ^iM)I[S-^\] ■■■]), V^(W) = n^7(n.)n^(^^-^.)' (1-3) 

ni,...,njvf i=l i=l i<^j 

with d{n) = sm{7m/N), g{n) oc exp{2TTi{J / N)n) and M = J = N/2 has been found to 
reproduce the ground state energy of the nearest neighbour Heisenberg chain of sites with 
remarkable accuracy^. Furthermore, ( |1.3|) captures the essence of the spin spin correlations 
in this system. States of the form (|1.3| ) span a large part of the Hilbert space including the 
ground state of the Haldane Shastry (HS) spin chain with long range exchange interactions 
Jki (x 1/ sin^(7r(A; — /)/A^)lll. The spectrum of this model gives a representation of the k = 1 
SU(2) WZW conformal field theor}^ (note that the next nearest neighbour interaction in 
the HS model is very close to Oc)- 

The spin Peierls system and ladder like models on the other hand have been studied 
using short range RVB states as ground states and soliton states interpolating between the 
two singlet bond configurations ( |1.2|) for the excitation^. 



In the following section we introduce variational states depending two free parameters 
and prove that they are eigenstates of the total spin of the system. In Section |11 B| we show 



that these states contain the valence bond states ( |1.2| ) and the Jastrow state (|1.3| ) as certain 



limits. Their relation to Gutzwiller projected states of spin-| fermions is discussed in [11 C| . 
Finally, we apply the variational states to the spin-Peierls system (|1.1| ) in Section |ITJ 



II. THE VARIATIONAL STATES 



A. Construction 



Let Ci and C2 be two parallels with polar angles 6 and it — 6 respectively on the 
two-dimensional sphere S*^. For even we choose N/2 equidistant points labeled by 
even numbers 0, 2, 4, ... , (A^ — 2) on Ci and similarly N/2 points labeled by odd numbers 
1,3,5,..., (A^ — 1) on C2 such that their azimuthal angles satisfy 

271 

— {l-2p)=^2i+i-V2i (2.1) 

for i G {0, . . . , (A^ — 2)/2}. This defines a lattice on S'^ characterized by the parameters 9 
and p consisting of two identical linear sublattices and obeying periodic boundary conditions 
in the azimuthal direction. The lattice points a can be described in terms of polar coordi- 
nates on the sphere, i.e. a = {cos (p a sin 9 a, simp a sin 9 a, cos 9 a). Alternatively, they may be 
parametrized by spinor components Ua = cos{9a/2) exp(z(y9a/2), Va = sin{9a/2) exp(— z(y9a/2). 
Given two points a,h E S"^ we construct the spinor product 

d{a, h) = {uaVh - UbVa) . (2.2) 

that has the following properties: 

\d{a,b)\' = \\\i<^-m\ 

... , _J sin((g, + g,)/2)sin((y;.-y.,)/2) ^^ 

argiala, 0)1 = tan ^ , , , z-rr \ ? (2.3) 

^ ' " \sin{{9a-9,)/2)cos{{^a-Vb)/'2)l ^ ' 



where 11 ■ 11 stands for the euclidean norm in R' 
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Now let each lattice point carry a spin S„ with S'„ = |. In the Hilbert space of these 
spins we consider the following states {d given by (p.2|)): 

i^™(^^,p))oc y: n^^^'"" fn^'K^.) 1 !■■■!)• (2.4) 



\{n,}\=M i 



V i<3 

Here | | ■ ■ ■ t) is the ferromagnetically ordered state with all spins parallel. The sum extends 
over all possibilities to select M out of lattice sites and invert their spins. The sets of the 
selected sites are denoted by {rii}. Each of those spin-product states contributes with an 
amplitude that consists of a product of single-site phasefactors and a product of Jastrow- 
like two-site factors, depending on the underlying lattice, namely the parameters 9 and p. 
The spinor components on the nth lattice site may be rewritten up to a common factor 
as Un oc exp(+i((7r/A^)(?T, + (— 1)"^)), f„ oc exp(— 2((7r/A^)(ra + (— with w = p + in 
and K = —{N/27!') ln(cot(^/2)). Hence, the states defined in (|2.4| ) depend analytically on a 
complex parameter expressing a dimerized structure of the lattice and will be referred to as 
dimerized Jastrow states (DJS). Periodic boundary conditions require the parameter J to 
be taken integer. Requiring the DJS to have a definite total spin restricts J even further: 
Clearly, the ( p.4|) are eigenstates of the z-component of the total spin with eigenvalue N/2 — 
M. To show that they are eigenstates of the total spin we rewrite as J2 ^n+J2n<m{Sn + 
+ 25^5^). Obviously E„ \i;^'^'J{9, p)) = fiV \i;^''^'^'J{9, p)) and 

n<m 

= ^ {M{M- 1) + {N- M){N- M- 1) - 2{N- M)M) \i)^^^^\9, p)). (2.5) 

For integer J with M < J < N — M this gives together with the results of appendix 0: 

\r^''''{9,p)) = (y - M)(| - M + 1) |^^'^'^(0,p)) . (2.6) 

Hence, under the above condition the state I'lp^''^^''^ {9, p)) is an SU(2) highest weight state 
with total spin 5* = N/2 — M for arbitrary values of the parameters 9, p. In particular, 
choosing M = J = N/2 one obtains a singlet, M = N/2 - 1, J G {N/2 - 1, N/2, N/2 + 1} 
gives three triplet states with = 1 etc. 



B. Limiting cases 



From (2^) one finds that the amphtudes of the wave functions (|2.4|) will be complex in 
general. These 'chiral' spin states are believed to arise in two dimensional spin systems due 
to frustrating interactionslll (see also Section |111 A| below). Choosing w G R, i.e. 6* = |, 



all lattice sites lie on the equator of the sphere and the Jastrow-factors in (|2.4| ) become all 
real: now the relative phases of the amplitudes are determined by the one-site factors alone. 
In this case, some of the states discussed in the introduction can be obtained by properly 
choosing p: 

(i) w = 0. Provided that M - 1 < J < iV - M + 1, the DJS 

l^^f'^) := l^™(f,0)) (2.7) 



are eigenstates of the Haldane-Shastry modelEl 



N-l 



„:tlo (-sm(f (n-m)))2 



with the eigenvalues 

EnM,J = ^(^)'(f - 1) + ^(^^' - 1) - 3MJ(iV - J)) (2.9) 

The ground state is given by the singlet, the first excitations by the triplet with M = 
N/2 - 1,J = N/2. 

In the thermodynamic limit the two spin correlation functions have been evaluated exactly 
by Gebhard and VollhardtS. For the singlet state of ( ^.7[ ) they have obtained 

hm I S.S.+„ I = (-1)"^^^ (2.10) 

N^oo 4 Tin 

This result is based on an alternative formulation of the DJS that is discussed below. 

(a) w = \. Here an even-numbered lattice site coincides with the next higher site. It is 

l^"--(^,^))oc|V^?],). (2.11) 
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The singlet of the DJS at w = 1/2 is a valence-bond-state: next neighbours are coupled to 



singlets. is easily proved by calculating the scalar product of both sides. Similarly, 



we have | ^"'^ (|, — |)) oc l^/'y^) 



I r^'^-'^^i^, t° ' ^] • ■ ■ I ^^''^ ■ • 2) ; {N- 1)] (2.12) 

DJS with higher total spin cannot be cast in a comparably simple form at w = 1/2 due to 
the more complicated structure of the remaining amplitudes. 

The two valence-bond singlets span the ground state space of the Majumdar-Ghosh modeli 

N-l ^ 

HmG = X! (^n^n+l + -S„S„_|_2) (2-13) 
n=0 ^ 

Recently, Nakano and Takahashi have generalized this hamiltonian to a variety of models 
with interactions of arbitrary range that have the same property^. In the thermodynamic 
limit these models have a finite gap for excitations over the twofold degenerate ground 
state leading to spin spin correlations that decay exponentially at large distances. This is in 
contrast to the Haldane-Shastry model which has no gap and according to ( p^.lOD correlations 
decaying algebraically. 

Hence the DJS with Q = -n jl and M = J = N/2 interpolate between a 'resonating valence 
bond' singlet for p = and the nearest neighbour 'valence bond solid' described by ( |1.2|) 
for p = |. Despite their essentially distinct properties these states have a remarkably large 
overlap. In Figure |l] the squares of the overlaps as computed numerically are plotted for 
systems up to = 24. At = 24 the singlet-subspace has the dimension 208012, so that 
a square overlap of 0.17 should be regarded as quite large. 

It would be interesting if one could find a hamiltonian interpolating between one of the 



, ■ — -1 
,vr 2 ^ 



models given in Ref. |20| and the Haldane-Shastry model, so that the singlet-DJS is always 
the exact ground state. 

(Hi) 6 ^ 0. For 6* — > and finite p the chains Ci,C2 are drawn to the poles of the sphere. 
In ( |2.4| ) all amplitudes and consequently the normalisation sum become zero. In this limit 
tensor products of Haldane-Shastry-typed states arise. In the simplest case, the singlet-DJS, 



one obtains for even N/2: 

N N H K H iL N N 

lim |^^'T.T(^^,p))oc \^P^s"')®\i^Hs"'), (2.14) 

o — 

The first factor in the tensor product refers to Ci, the second to C2. Other DJS can also be 
examined by exphciting their dependence on the parameters 6 and p. 



C. Relation to the Gutzwiller wave function 

It is worthwhile noticing that the DJS can also be formulated analogous to the Gutzwiller 
wave functionill. 
Let 

ctna^ ^na be cauouical creation and annihilation operators on the nth site for spin- 2- 
particles with S^ = a and | ) the vacuum of the system . 

For e.g. M = J = N/2 one can use the following construction: With the definition 

1 N-l 

btH = E e^^("+(-^)"-)a+^ (2.15) 



n=0 

we have 

|^^'T.f(^,p))ocP6+^H---6+^ ^{w)bl^{w)...bt^ ^H|0) (2.16) 

with kn = + § + andV being the Gutzwiller projector excluding double occupancies, 
i.e V = n^=o^(l~'^nT^nj,)5 rina = a^^a„o-. Eq. ( |2.16|) can be proved using similar arguments as 



in Refs. |T5|. So the DJS are deformed Gutzwiller wave functions characterized by a complex 
parameter w. Analytical results concerning the two-site correlations have been found for 
w = dliii. For non-zero w the dimerization leads to complications that cannot be resolved 
following these methods. 

An obvious generalisation is to construct wave functions of the above type with arbitrary 
filling. Starting from states without double occupancies on the equidistant lattice as con- 



sidered in Ref. ^ one can introduce double occupancies and dimerization while keeping the 
Jastrow form of the states^. Again they can be constructed to be heigest weight states of 
the total spin. Since empty and doubly occupied lattice sites do not contribute to the spin, 
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there are large additional possibilities of varying the wave functions in the case of general 
filling and strength of the Gutzwiller projection. 



III. APPLICATION TO THE SPIN-PEIERLS SYSTEM 



We now use the DJS as variational ansatz for the low lying states of the model ( |1.1| ). As 



mentioned above the analytical methods of Refs. p!9| , |22| cannot be applied to the dimerized 
system, hence the results presented below were obtained by numerical evaluation of the 
relevant matrix elements for system sizes up to = 26 lattice sites. 

A. Ground state properties 

For the ground state of the model we have used the singlet-DJS (M = J = N/2) as a 
variational ansatz. The numerical results can be summarized as follows: 

For 2a + 6 < 1 the expectation value of the hamiltonian is minimized by real w, i.e. 
6 = 7t/2 with p varying from to 1/2, 2a + 6 = 1 corresponds to a valence bond state 
(w = 1/2)0. On the other hand, for 2a + 5 > 1 we find p = 1/2 and varying 6. For 
a = S = the minimum is given by the Haldane-Shastry ground state {w = 0). The exact 
ground state energy per spin is known to be — In 2 + 1/4 ~ —0.443147 for an infinite system, 
the variational value is — (3/4)Si(7r)/(7r) ~ —0.442177 from (|2.10|) . In Table | we present the 



extrapolation to an infinite system of the variational ground state energies for various values 
of a and 6. Comparison with data from numerical diagonalisation obtained by Chitra et 
al.ii shows excellent agreement: the variational energies per spin differ by only about 10~^ 
from diagonalisation values except for the last row of table |. In general, the quality of the 
singlet-DJS as a variational ansatz is found to decrease for values of a exceeding 0.5. Hence, 
we have found an effectively one-parametric variational wave function that is an excellent 
approximation of the true ground state of the model considered within a large area of its 
parameters, in particular for any antiferromagnetic nnn couplings a with 2a + 6 < 1. Note 



that in this case the optimum singlet DJS (|2.4| ) has real amplitudes. For 2a + 6 > 1 they 



become properly complex: this case therefore is called a chiral phase. 

For 6 the translational invariance of the system is explicitly broken leading to a dimerized 
variational ground state with finite p. 

For 6 = and small values of the nnn coupling < a < a* the variational energy is 
minimized for w = 0. A transition to a dimerized ground state is observed at the 'critical' 
value of a* ~ 0.2716 ± 0.0002 for — > oo. For a > a* p becomes finite. In absence of 
the alternating term the Hamiltonian {6 = 0) the system is invariant under translation by a 
single lattice site, while the states ( p.4|) do not have such a symmetry for w 7^ 0. This leads 
to consider the ansatz 

^r / N , : AT N N ,71 ^ ^ , N N ,7r 

|^^(p))0c|^^'T-(_,p)) +|^A^>-.-(_,_p)) (3.1) 

with N chosen even. This construction corresponds to the lattice momentum of the ground 
state of finite systems obtained in Ref. ^ With the states ( p.l| ) one observes slight, but 
numerically significant improvements of the ground state energy as shown in table |l|. Fur- 
thermore, the parameter p remains zero for < a < a** with a** ~ 0.1737 ± 0.0002, which 
is much smaller than the a* mentioned above. In Figure ^ we present the ground state 
energy per spin as a function of a, in Figures |^,^ nn and nnn correlations calculated within 
the states (|3.1| ) are plotted. These diagrams agree very well with the corresponding figures 
given in Ref. ^ 

Spin-spin correlations beyond those entering the expression for the ground state energy 
show the correct long distance asymptotics as long as p = (|2.10|) , which corresponds to the 
massless regime a ^ ^, 6 = described by a conformal field theory. We find that non-zero 
p lead to a suppression of long-range correlations. The system sizes that we have analyzised 
numerically do not allow, however, to study the dependence of the correlation length on p. 

B. Excitations 

The results of the previous section suggest to use DJS with higher spin as variational 
ansatz for excitations of our model. Unfortunately the situation is not as clear as before. 
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Here we concentrate on the case S = and mention that eigenstates of the Haldane-Shastry 
model provide a good description of low-lying states for < a < 0.3. 

From analytical and numerical studies this model is known to be gapless for a < ac- As 
predicted by conformal field theory, at the 'conformal point' ac is defined by the occurence 
of many degeneracies in addition to the usual SU(2) symmetry. Okamoto and Nomura 
examined the spectrum of (|1.1|) for 5 = in finite systems by numerical diagonalisationi. 
They found a linear dependence of the energy of the ground state and the two first ex- 
citations (triplet, singlet) on a in the above intervall. From the condition, that the two 
excitations should degenerate, they obtained a precise value for the conformal point as 
etc = 0.2411 ±0.0001. Within the concept of DJS, there is only one singlet corresponding to 
the groundstate, but in the HS case {w = 0) an additional singlet can be derived from the 
triplet excitation by using the Yangian symmetry of this modelil. 

The Yangian of the Haldane-Shastry model is generated by the total SU(2)-spin or 'level- 
operators' Qq = I] 5" and the 'level-l-operators' Q° = X^m^n cot(;^(m — n))e°'^^ S^S^. 
We do not repeat any details of this symmetry algebra and its representations here. In 
the HS model the lowest excitation for an even number of spins is given by the triplet- 
DJS \'^h's^'^~^'^^'^)- This state degenerates with a singlet that can be obtained by applying 
Qi = Qi ~ '^Q\ oil \'4'h's^'^ g^j-^^ projecting onto the singlet space, i.e. 

N i^ — 1 H- 

|^f5)oc(2-Q2)Qrl^H5^''^). (3.2) 

As mentioned before, the minimum expectation value of the hamiltonian on the line 5 = 
is given by the HS ground state for a < a* . So we evaluated the hamiltonian also for the 
two excited HS states above. Note that in this ansatz there is no variational parameter 
included, because ( p. 21 ) is not defined for w 7^ 0. In Figure || the difference of these energies 
with the ground state value given in Figure |^ are plotted. For a < 0.3 there is a good 
agreement with the corresponding data of Ref. ^ At some a = a the singlet and triplet 
energies coincide. Table |T| shows values for a in finite systems, which can be extrapolated 
to ^ 00 smoothly giving a = 0.2368 ± 0.0002. This is remarkably close to the value 
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= 0.2411 ± 0.0001 obtained from numerical diagonalisation. 



IV. CONCLUSIONS 

We have constructed a family of variational states which contain excellent approximations 
to the ground state of the model ( |1 . Ij ) for a large range of its parameters. The states are 
strictly spin-singlets and can be characterized by a single complex parameter. Moreover, in 
Section |TTl| we have demonstrated that for 6 = and < a < 0.3 groundstate and lowest 
excitations of our model can be desribed amazingly well by eigenstates of the HS model. 
For 6^0 the variational ansatz reproduces the explicit dimerisation of the model.0 For 
vanishing bond alternation 6 = the ground state of the finite system has a well-defined 
lattice momentum kg G {0, 7r}@ while the lowest excitations (triplet, singlet) have kt, kg* = 
TT — kg (mod 27r). In the thermodynamic limit the excitation gap of the triplet vanishes for 
< ttc; while for larger a the singlet degenerates with the ground state of the finite system. 
The momenta of these two states differ by tt. Consequently, in this two dimensional ground 
state space linear combinations exist that represent dimerized states. This is expressed 
within our variational approach by the fact that even for finite systems dimerized states 
(with p 7^ 0) are good approximations of the ground state. The corresponding translational 
invariant states ( p.l|) lead to small improvements to the ground state energy, but reproduce 
the spin correlations between nearest and next nearest neighbours very well. 

As mentioned above a generalization of this variational states to similar systems away 
from half filling is straightforward. This may allow for similar studies of the Zn-doped 
compound Cui„a.Zn^Ge03ii when the Zn-sites are treated as static spin-0 sites in the spin 
chain. 
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APPENDIX A: OFF-DIAGONAL MATRIX ELEMENTS OF THE DJS 



To compute the action of J2n<m{^n^m + ^m) on the DJS (p^ ) we define 

lx(W)) = n^nJT---T) (Ai) 

i 

with \{ni}\ = M. Introducing z = e*^ and = i^z^'^)^^^^'^' we find 

{x{{n^})\r^''^\e,p)) 



E-'^"E( n 

n=l i iT^j 



II lyTL'i lyTliA 1 — TX i ■y Th -j 

i^j H-n.Z Z J Z Z J 



^ 1 J_^ni+2n+l _ n 



_|_ ^J(2n+1) / TT H-^iZ Z 

|S2N 

^ni+2n+l ^n,- ^— n^ — 2n— 1 .v— 

^ ) A2 

The two sums originate from hopping terms between spins separated by an even and odd 
number of lattice sites respectively. Combining these terms ( |A2|) can be rewritten as 

n<m 

where P is a polynomial in its arguments with coefficients independent on n. Only powers 
(±2;)"^ with — (M — 1) < A; < (M — 1) arise in the expansion so that for integer J 



with M < J < N — M the sum taken over each term of P vanishes leading to 

E (S^S- + S~S^) p)) = -M p)) , (A4) 



n<in 
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TABLES 



TABLE L 


Variational ground state enerj 


^ies per spin for N — 


> oo. The last column shows 


numerical dat 


a from Ref. p^. 






Q 


6 






0.2411 


0.00 


-0.40160 


-0.401866 


0.2411 


0.04 


-0.40793 


-0.409051 


0.2411 


0.16 


-0.44237 


-0.442862 


0.2411 


0.32 


-0.49673 


-0.496844 


0.2500 


0.00 


-0.40008 


-0.40045 


0.2500 


0.35 


-0.50721 


-0.50727 


0.2500 


0.20 


-0.46242 


-0.46329 


0.4000 


0.50 


-0.56550 


-0.56611 


0.4800 


0.10 


-0.41272 


-0.41281 


0.5500 


0.10 


-0.41519 


-0.41610 


0.2000 


0.80 


-0.67601 


-0.67613 


0.3000 


0.80 


-0.67896 


-0.67966 


0.4800 


0.80 


-0.68897 


-0.69256 



TABLE n. Ground state energies per spin at (5 = and oo. The second column shows 

variational results obtained from ( |2.4[ ), the third from (3.1). The data in the last column is taken 
again from Ref. p5l 



a 



i?g-(a,0,oo) 



T-^diag 



/N 



0.20 
0.25 
0.30 
0.40 



-0.40852 
-0.40008 
-0.39158 
-0.37919 



-0.40880 
-0.40045 
-0.39240 
-0.37924 



-0.40885 
-0.40045 
-0.39284 
-0.38028 
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TABLE III. Location of the degeneracy a between singlet and triplet excitation for different N 

N a 

10 0.23967 

12 0.23876 

14 0.23820 

16 0.23785 

18 0.23757 

20 0.23739 

22 0.23722 

oo 0.2368 
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FIGURES 

FIG. 1. Overlap 0)|V'^'^/2'^/2(|, i))|2 between the ground state of the Hal- 

dane-Shastry model and the valence-bond-state. 
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Figure 1 from H. Frahm and J. Schliemann 
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FIG. 2. Variational ground state energy per spin at S 



= as a function of a. 




21 



FIG. 3. Nearest neighbour spin correlations {S^S^_^_l) = |(S„S„-t-i) in the variational ground 
state of (|l.l|) for 5 = as functions of a calculated from (|3.lD for different system sizes N. 
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FIG. 4. Next nearest neignbour spin correlations {SnSn+2) — \{^n^n+2) of the variational 



ground state of (|l.lD for 5 = as functions of a calculated from (|3.l|) for different system sizes A'". 
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FIG. 5. Energy gap between the variational ground state and the lowest triplet Agt and the 
singlet (|3.2| ) A^^* for (5 = as a function of a for N = 20. The degeneracy determines our estimate 
of the conformal point ac- 




